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SUMMARY 

The  flow  of  fluids  with  and  without  heat  sources  and  subject  to 
body  forces  between  two  plane  parallel  surfaces  which  are  oriented  in 
the  direction  of  the  generating  body  force  is  analyzed  under  the  condi- 
tion that  the  temperature  vary  linearly  along  these  surfaces.  It  is 
found  that  a modified  Rayleigh  number  (product  of  the  reciprocal  of  the 
ratio  of  specific  heats  and  the  Prandtl  number  Pr  and  the  modified 

Pfyd 

Grashof  number  Gr/J  as  well  as  a parameter  = Pr  Gr^  — — is  of 

CP 

significance  in  this  problem;  where  B is  the  volumetric  expansion 
coefficient,  fx  is  the  negative  of  the  X-component  of  body  force  per 
unit  mass,  d is  the  characteristic  length,  and  Cp  is  the  specific 
heat  at  constant  pressure.  Solutions  of  this  problem  are  obtained  in 
terms  of  "universal”  functions  which  are  tabulated  for  simple  applica- 
tion to  specific  cases.  Representative  velocity  and  temperature  dis- 
tributions from  which  detailed  study  of  the  heat  transfer  is  made  are 
then  computed.  When  the  ratio  of  CK^  (where  C is  related  to  the  mass 
flow)  to  the  Rayleigh  number  is  of  unit  order  of  magnitude,  the  effects 
of  aerodynamic  or  frictional  heating  can  be  appreciable.  Asymptotic 
solutions  (for  large  values  of  the  Rayleigh  number)  which  render  the 
computations  simple  are  also  presented. 

Comparison  of  the  results  from  the  method  given  herein  with  those 
obtained  elsewhere  in  an  approximate  manner  for  a special  case  simu- 
lating the  natural-convection  flow  of  fluids  with  heat  sources  in  a com- 
pletely enclosed  region  shows  that  the  approximate  method  is  suffi- 
ciently accurate  for  problems  in  which  the  modified  Rayleigh  number  is 
less  than  104. 
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INTRODUCTION 

In  recent  years  the  transfer  of  heat  to  and  from  enclosed  or  par- 
tially enclosed  regions  by  means  of  natural  convection  or  by  a combina- 
tion of  natural  and  forced  convection  has  taken  on  new  significance  in 
the  fields  of  aeronautics,  atomic  power,  electronics,  and  chemical  en- 
gineering. Most  of  the  information  on  these  modes  of  heat  transfer  under 
such  conditions  is  of  a semiempirical  or  specialized  nature;  relatively 
little  detailed  information  exists  for  internal  natural-convection  flows. 

In  reference  1 there  appears  one  of  the  few  attempts  to  determine  theo- 
retically the  velocity  and  temperature  distributions  in  detail  and  hence 
the  heat  transfer  for  an  internal  flow  problem  of  thiB  kind.  In  that 
reference  a solution  was  found  for  the  fully  developed  flow  of  fluids 
with  and  without  heat  sources  between  two  long  parallel  plates  with  con- 
stant wall  temperatures  (where  one  could  be  different  from  the  other) 
oriented  in  the  direction  of  the  generating  body  force.  The  information 
obtained  therein  is  of  practical  value  in  connection  with  fully  developed 
flows  subject  to  body  forces  where  the  surfaces  are  maintained  at  uniform 
temperatures.  Lighthill  (ref.  2)  employs  integral  methods  to  study  the 
natural-convection  flow  in  tubes  with  either  one  end  or  both  ends  closed 
and  with  constant  wall  temperatures,  and  in  reference  3 approximate 
superimposed  free-  and  forced-convection  flows  are  obtained  for  short 
channels  and  pipes. 

As  the  next  step  in  the  study  of  natural  convection  or  combined 
natural  and  forced  flows  in  confined  spaces,  consideration  is  here  given 
to  the  configuration  of  reference  1 with  the  exception  that  the  thermal 
boundary  condition  specified  is  that  the  surface  temperatures  vary 
linearly  along  the  plates  or  surfaces.  (One  surface,  however,  may  be  at 
a different  local  temperature  from  the  other  but  the  slopes  of  the  tem- 
perature distributions  on  each  surface  are  taken  to  be  equal.)  The 
analogous  forced-convection  problem  is  treated  in  references  4,  5,  and  6. 
The  present  problem  simulates  several  important  physical  occurrences  of 
this  phenomenon;  for  example,  it  could  represent  the  case  where  the  out- 
side of  the  channel  formed  by  the  plates  is  cooled  (or  heated)  by  a 
counterflow.  In  addition,  the  present  problem  represents  a more  general 
case  than  was  considered  in  reference  1,  since  here  the  temperature  will 
no  longei  be  restricted  to  be  a function  of  the  transverse  coordinate  alone 
and,  hence,  energy  convective  as  well  as  mass  convective  effects  will  be 
included. 

The  solution  is  obtained  in  terms  of  functions  which  depend  on  only 
one  of  the  several  associated  dimensionless  parameters,  and  these  func- 
tions are  tabulated  so  that  specific  cases  can  be  easily  computed.  Solu- 
tions for  pure  natural  convection  and  for  superimposed  natural  and  forced 
convection  are  shown  to  be  essentially  identical.  Representative  velocity 
and  temperature  distributions  are  also  presented,  and  the  effects  of 
frictional  or  aerodynamic  heating  on  the  flow  and  heat  transfer  are 
discussed. 
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A special  case  which  simulates  a completely  enclosed  region  In  which 
there  Is  no  net  mass  flow,  the  walls  are  at  the  same  temperature,  and 
heat  is  generated  uniformly  "by  heat  sources  is  treated,  and  detailed 
velocity  and  temperature  profiles  are  obtained  from  which  the  heat  trans- 
fer is  determined.  This  special  case  was  treated  in  an  approximate 
manner  in  reference  7. 

Consideration  is  also  given  to  the  problem  of  convective  inversion, 
that  is,  to  the  cases  where  the  modified  Grashof  number  changes  sign. 

It  is  shown  that  convective  inversion  due  to  changes  in  the  body  force 
direction,  to  changes  in  the  sign  of  the  volumetric  expansion  coeffi- 
cient, or  to  changes  in  the  sign  of  the  axial  temperature  gradient  alters 
the  character  of  the  problem,  because,  it  is  believed,  under  these  con- 
ditions the  flow  becomes  unstable  because  of  heating  from  below  (see 
pp.  104  to  107,  ref.  8). 


ANALYSIS 

Formulation  of  the  Problem 

The  study  to  be  made  here  is  that  of  the  laminar  fully  developed 
flow  of  fluids  with  and  without  heat  sources  and  subject  to  a body  force 
between  two  plane  parallel  surfaces  open  at  both  ends  and  oriented  in 
the  direction  of  the  generating  body  force  (see  fig.  l).  It  is  further 
v specified  that  there  shall  be  linear  (with  equal  slopes)  temperature 

variations  along  the  walls  but  that  the  walls  need  not  necessarily  be  at 
the  same  temperature.  The  flow  is  assumed  to  be  parallel  to  the  axis  of 
! ’ the  channel  (that  is,  the  only  nonvanishing  velocity  component  is  the  one 

in  the  longitudinal  direction)  and  in  addition  it  is  assumed  that  the 
i ■ physical  properties  (for  example,  Cp  and  |i)  of  the  fluids  are  constants 

and  that  the  essential  influence  of  the  density  changes  on  the  flow  is 
taken  into  account  by  the  introduction  of  the  volumetric  expansion  co- 
efficient in  the  body  force  term  (that  is,  the  other  influences  of  vari- 
able density  and  the  variation  of  the  expansion  coefficient  with  tempera- 
ture are  negligible).  Discussions  of  the  Justification  of  the  assump- 
tions can  be  found  in  references  1,  9,  and  10. 


Under  the  conditions  stated,  the  basic  equations  with  body  forces 
included  expressing  the  conservation  of  mass,  momentum,  and  energy  (see 
ref.  l)  become,  respectively, 
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dP 

3y 


0 


(3) 


^2T»  d^T*  ^vTIdT*  £ /au\2  Q 

ax2  ay2  " k sr - k \5f ) " e 


(4) 


(See  appendix  A for  a complete  list  of  the  symbols  used  herein.)  With 
cognizance  taken  of  equations  (l)  and  (3)  the  above  system  reduces  to 


d^  1 ( . dp\ 

dY2  " M rX  «/ 

a2!*  a%*  ^vti^t*  j± /«j\8  q 
ax2  ay2  " k sr  " k \&)  “ e 


(5) 

(6) 


where  the  velocity  is  a function  of  the  transverse  coordinate  Y only 
and  the  pressure  is  a function  of  only  the  longitudinal  coordinate. 

The  body  force  term  in  equation  (5)  can  be  written  as  a buoyancy 
term  by  introducing  the  volumetric  expansion  coefficient  f)  in  a manner 
similar  to  that  described  in  appendix  B of  reference  1.  Equation  (5) 
then  becomes 


§ ♦ £ (T«  - T.0.)  . i (f  ♦ P^*)  (7, 

where  the  subscript  wq  refers  to  the  surface  at  Y * 0 (see  fig.  l). 

The  boundary  conditions  associated  with  this  problem  are  that  the 
velocity  at  the  walls  must  vanish  (the  no-slip  condition  for  viscous 
fluids)  and  that  the  temperature  must  vary  linearly  along  the  walls. 
(Note  that  the  latter  condition  implies  that  the  temperature  gradients 
along  the  walls  and  hence  the  axial  heat  flux  along  the  walls  must  be 
constant.)  In  order  to  satisfy  the  temperature  conditions  and  equation 
(7),  the  temperature  must  be  of  the  form 

T*(X,Y)  - AX  + T(Y)  (8) 

Mathematically,  the  boundary  conditions  are  formulated  as 


U(0)  - U(d)  - 0 
T*(X,0)  - AX  + T(0) 


(9) 

(10) 
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T*(X,d)  - AX  + T(d)  (11) 

Substituting  equation  (8)  into  equations  (7)  and  (6),  respectively, 
yields 


and 


where 


(12) 


(13) 


0 - T(Y)  - T(0) 

Since  the  left  side  of  equation  (12)  is  a function  of  Y alone 
and  the  right  side  is  a function  of  X alone,  it  is  clear  that  each 
side  must  be  equal  to  a constant.  Thus,  equation  (12)  can  be  written  as 


dY2 


(14a) 


or 


5 (§  + V*)  - c C«») 

The  temperature  now  appears  in  equations  (13)  and  (14)  as  0 which 
is  independent  of  the  longitudinal  coordinate  X,  and  hence  the  thermal 
boundary  conditions  can  be  written  as 

0(0)  - 0 

0(d)  =»  T(d)  - T(0)  E 0Wl  (15) 

To  nondimen s iona li z e equations  (14)  and  (13),  let 
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ku 


pefxd 

Y - yd 


g m ti SI T 

0 p2p2fx2d4  " KA 


(16) 


where  KA  * PrGrA 


/PfXd\ 


is  the  new  dimensionless  parameter  presented 

\ °p  / 

in  reference  1,  but  here  the  Grashof  number  is  modified  because  it  is 
based  on  Ad  rather  than  on  a temperature  difference,  that  is, 
pfvd4A 

GrA  - -jr-- 


Hence,  equations  (14),  (13),  (9),  and  (15)  become 


u"  + T « CKA 

(17) 

- Ra  u + (u*)2  + aXA  a 0 

(18) 

u(0)  - u(l)  = 0 

(19) 

T (o)  - 0 

(20) 

KA*Wi 

Td)  ■ *, 

(21) 

where  the  primes  denote  differentiation  with  respect  to  y, 

Ra  = (l/r)  PrGrA  is -a  modified  Rayleigh  criterion  or  number  (see  p.  105, 

ref.  8),  a ■ Qd/kA  is  the  dimensionless  heat  source  parameter,  and 
C *»  C /^Prd3/cpAKA.  Eliminating  T between  equations  (17 ) and  (18) 
results  in 

uiv  - (u')2  + Ra  u - aKfr  ■ 0 (22) 
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with  the  boundary  conditions 

u(0)  ■ u(l)  = 0 

(23) 

^ s 

CO 

II 

a 

& 

(24) 

u" (l)  = (c  - eW]/Ad)  Ka  = mCKA 

(25) 

where  m «■  1 - 0w/CAd.  The  constant  C which  appears  as  a parameter  in 
the  boundary -value  problem  described  by  equations  (22)  to  (25)  merely 
specifies  the  temperature  level  (see  eq.  (17))  of  the  problem.  In  order 
to  define  completely  the  temperatures  and  velocities,  this  constant  must 
be  related  in  some  way  to  the  physics  of  the  problem.  From  equation  (14b) 
it  can  be  seen  that  C could  be  determined  from  the  pressure  gradient 
along  the  channel;  that  is,  C is  essentially  connected  with  the  end 
conditions  to  which  the  channel  is  subject.  Since  the  pressure  gradient 
may  not  be  known  a priori,  in  the  subsequent  section  dealing  with  the 
solution  of  the  present  problem  C will  be  related  to  the  end  conditions 
by  the  mass  flow  in  the  channel,  which  remains  invarient  over  the  entire 
length  of  the  channel. 

Note  that  solution  of  the  preceding  boundary-value  problem  will 
yield  velocity  and  temperature  distributions  for  both  natural-convection 
and  combined  natural-  and  forced-convection  flows.  The  forced- 
convection  pressure  gradient  merely  alters  the  magnitude  of  the  constant 
C.  A discussion  of  such  a superimposed  flow  problem  under  special  con- 
ditions is  given  in  reference  11. 

Several  interesting  observations  can  now  be  made  concerning  these 
equations.  First,  comparison  of  equation  (22)  with  the  corresponding 
equation  in  reference  1 shows  that  they  are  identical  except  for  the 
third  term  in  equation  (22),  which  does  not  appear  at  all  in  the  equa- 
tion in  reference  1.  This  term  stems  from  the  convection  term  in  the 
energy  equation;  hence,  in  the  present  problem  energy  convection  effects 
will  be  included.  The  energy  convection  term  vanished  identically  in 
reference  1 because  of  the  assumption  that  the  velocity  and  temperature 
profiles  were  independent  of  the  axial  coordinate.  Second,  since  the 
convection  term  appears  with  a coefficient,  another  dimensionless  param- 
eter (Ra,  the  modified  Rayleigh  number)  is  associated  with  this  problem 
and  its  influence  on  the  results  must  be  studied.  Finally,  the  condi- 
tions of  the  problem  require  that  the  temperature  be  of  the  form 
T*  » AX  + T(y)  (see  eq.  (8)),  and  hence  the  longitudinal  heat  flux  is 
everywhere  constant. 
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The  boundary-value  problem  stated  in  equations  (22)  to  (25)  can  be 
written  in  more  convenient  form  by  defining 

7 ■ Ra  u - aKA  (26) 

Hence 

v*7  + Ra  v - ^ (v1  )2  ■ 0 (27) 

7(0)  - 7(1)  - - oXA  (28) 

7"(0)  = + Ra  CKa  7"(l)  - + m Ra  CKa  (29) 


Solutions  of  the  Boundary-Value  Problem 

Equation  (27)  is  nonlinear  (the  nonlinear  term  is  due  to  the  fric- 
tional or  aerodynamic  heating)  and  therefore,  as  in  reference  1,  a method 
of  successive  approximations  will  be  employed  to  find  its  solution.  To 
this  end,  equation  (27)  is  written 

^ + Ra  ^ (7A-1>2  “ 0 (30) 


where  n - 0,  1 denotes  the  particular  term  in  the  approximation 
v ■ Vq  + v^  and  v*^  s 0.  Let 
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Equations  (30),  (28),  and  (29)  become,  respectively, 


viv  + Rv  - -(v*  , )2  . o 
n n r'  n-1' 

(32) 

voC_1)  " vo0-)  " * 

(33) 

V^(-1)  - vg(l)  - mJ 

(34a) 

Vn(-1)  » Vn(l)  - v”(-l)  - v£(l)  a 0 n + 0 

(34b) 

where  the  subscripts  now  denote  differentiation  with  respect  to  tj. 
Zeroth-order  approximation.  - In  the  zeroth-order  approximation, 

the 

not 

nonlinear  term  which  is  associated  with  the  frictional  heating  does 
appear  and  the  problem  then  consists  of  solving  the  equation 

vJV  + Rvq  a 0 

(35) 

subject  to  the  boundary  conditions 

v0(-i)  = vo(1)  3 X 

(36) 

v^(-l)  - JyfR  vj(l)  = mJ^R 

(37) 

For  simplicity  of  computations  the  solution  can  be  obtained  in  terms  of 
symmetric  and  antisymmetric  functions  of  tj,  depending  only  on  the  param- 
eter R (or  Ra),  by  setting 


% , J(m  - l)  , J(m  + l) 

v0  “ *v00  + — 2 — L V01  + — 2 — L V°2  (38) 

where  the  boundary  conditions  to  be  satisfied  by  the  various  vQj  (where 
j a 0,  1,  2)  are  given  in  the  following  table: 


Subscript 

v(-l) 

i 

v(l) 

v"(-l) 

v"(l) 

00 

1 

1 

0 

0 

01 

o 

0 

-*{r 

02 

0 

0 

VF 

Vr 
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The  functions  must  each  satisfy  equation  (35),  which  has  the  general 

solution  for  positive  values  of  the  Rayleigh  number 

vQ  = ax  exp(erij)  + &2  exp(-err))  + a3  e5q?(grt))  + a4  exp(-grrj)  (39) 

where 


e ■ i + 1 
g - i - 1 


and 

r = R1/4/^ 

Using  the  boundary  conditions  to  evaluate  the  constants  in  equation  (39) 
and  expressing  the  solution  in  terms  of  real  products  of  circular  and 
hyperbolic  functions  show  that 

v^  ■ 5 i = (cosh  r cos  r coBh  rq  cos  rTj 

w cosher  + cos^r  - 1 

+ sinh  r sin  r sinh  rq  sin  rq)  (40) 


vQ,  ■ 5 2~  (cosh  r sin  r sinh  rij  cos  rij 

cos  r - cosh  r 

- sinh  r cos  r cosh  ri)  sin  rn) 


(41) 


Vqo  = 2 2 (cosh  r cos  r sinh  rrj  sin  rr) 

v cosh  r + cos  r - 1 

- sinh  r sin  r cosh  rrj  cos  rn)  (42) 

Note  that  Vqq  and  Vq2  are  symmetric  functions  and  Vq^  is  an  anti- 
symmetric function  of  tj.  These  "universal"  functions  v^j  are  given 


for  various  values  of  Ra 


table  I. 
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First -order  approximation.  - In  order  to  include  the  effects  of 
frictional  heating  to  a first  order,  let  n ■ 1 in  equation  (32).  Thus 


(43) 


v£v  + Rvx  = | (^)2  - f(n) 

where,  in  view  of  equation  (38),  f (iq)  is  given  by 

tin)  - | [*2v<D02  + *J(m  + l)v00  v02  + i ^ 

+ J2(m  I).2  v.  2 + AJ(m  - l)v'  V + j2(m2  ~ 1)  V v'l 
4 02  A'  ■L;v00  01  2 01  02J 


'01 


(44) 


The  boundary  conditions  on  v^  are 

v1(-l)  = vx(l)  - v£(-l)  - v£(l)  = 0 
Once  again,  to  obtain  "universal-type"  functions  let 

V1  " x2yl0  + X' J(®  + 1^vll  + J '^m  4'  1" 


(45) 


'12 


J2(m  + l)2 


%,/  . J2(m2  - 1) 

M(m  - l)v14  + — i 1 


v13  + 


'14  T 2 15 


(46) 


In  view  of  the  form  of  equations  (43)  to  (46),  the  functions  v-^q,  v-^, 
v^g>  and  vi3  are  symmetric j and  v^4  and  v-^g  are  antisymmetric 
functions.  Each  of  the  functions  v-^  (k  = 0,  1,  ...  5)  will  be  a 

solution  of  equation  (43),  but  where  the  nonhomogeneous  term  is  only 
the  related  part  of  equation  (44)  (for  example,  the  first  term  on  the 
right-hand  side  of  equation  (44)  is  associated  with  the  v^q  solution 

and  so  forth). 


A particular  solution  of  equation  (43)  can  be  constructed  from 


(Vlk)p  mJ^  fk(  W °(Tl  " * 


■t 

* * 

*sr 


12 


NACA  TN  3141 


where  G Is  the  Green's  function  which  satisfies  the  homogeneous  part  of 
equation  (43)  and  the  conditions  G(0)  ■ G*  (0)  ■ G"(0)  ■ 0,  G"'(0)  ■ 1. 
These  arbitrary  conditions  were  applied  to  yield  a simple  form  for  the 
Green's  function.  Hence 


(v  ) 
v lk  P 


sinh  g r(rj  - {)  - g sinh  e r(ri 


Note  that  if  the  f^  is  a symmetric  function,  (vlk)p  i8  alBO  symmetric j 

and  if  fk  is  antisymmetric,  so  is  (vik)p*  From  the  boundary  conditions, 
equation  (45),  it  can  be  seen  that  V]_k  must  be  either  symmetric  or  anti- 
symmetric. Therefore,  if  fk  is  symmetric,  the  complementary  solution 
used  with  the  particular  solution  must  be  symmetric } and  if  fk  is  anti- 
symmetric, the  complementary  solution  must  also  be  antisymmetric.  These 
complementary  solutions  are,  respectively, 

(v  ) = E C06h  ?.rTl  + F cosh-SEH, 

' lk'cs  cosh  er  cosh  gr 


(v  ) s e Slnh  erT]  + F 
' lk'ca  sinh  er  sinh  gr 


(49) 


E 


/ [sinher  (1-^3 

r°o 


(50) 


F 


(51) 


The  v-y^  solutions  can  be  written  explicitly  in  real  form  (by  proper 

combinations  of  particular  and  complementary  solutions)  and  are  given  in 
equations  (Bl)  to  (B6)  in  appendix  B.  Values  of  v^  for  several  R 

(or  Ra)  are  presented  in  table  II. 


Velocity  and  Temperature  Distributions 

Now  that  the  various  v^k  are  known  explicitly,  they  can  be  inserted 
into  equation  (46)  to  yield  Vj_.  The  sum  of  Vq  (as  given  by  eq.  (38)) 
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and  v^  then  forms  the  solution  which  includes  frictional  heating  effects 

to  a first  order.  By  means  of  the  various  transformations  made,  the  solu- 
tions of  the  given  boundary-value  problem  (eqs.  (22)  to  (25))  or  the 
dimensionless  velocity  distributions  to  the  zeroth  and  first  approximation 
are,  respectively, 


"0 


KA  " 

Ra  -a 


v00 


(m  - 1)  Ra1/^  . (m  + l)  Ra1/2: 

+ 1 i v01  + * 1 


'02 


+ a 


(52) 


ui  - k 6‘ 


64  (vq  + V]_)  + aK 


(m  - l)2RaC2  . (m  + lJ^RaC5*  , „ l/2„ 

+ 1 f V12  + * 4 v13  - C«  - !)  oRa  / Gvu 


Uq  + 
s2n„«2 


ka2  r2 

64te  r v10~ 


aRa  1/2C(m  + 1) 


'11 


(m2  - l)RaC2 


-»] 


(53) 


where  the  vn.  (j  « 0,  1,  2)  are  given  in  equations  (40 ) to  (42)  and  the 


'lk 


OJ 

(k  - 0,  1,  2,  3,  4,  5)  are  given  in  appendix  B. 


In  principle,  higher-order  approximations  could  be  obtained  by  con- 
tinuing the  procedure  described.  However,  the  results  become  very  un- 
wieldy, Therefore,  beyond  the  range  of  applicability  of  the  zeroth-  and 
first-order  approximations  (that  is,  in  the  range  of  large  frictional 
heating  effects),  the  complete  boundary-value  problem  should  be  solved 
numerically^  some  discussion  of  these  numerical  results  relative  to  the 
zeroth-  and  first-order  approximations  will  be  presented  subsequently. 

To  determine  the  temperature  distributions,  recall  that  they  are 
related  to  the  velocity  distributions  by  equation  (17)  and  that 


v = 
00 


Ra 


V2 


v02  and  v02 


Ra 


1/2 


Vqq  so  that 


CK, 


n " 
"0 


Ka 


(C  ‘ [r02  + 2(m  ’ 1)CV°1  + (m  +21>C  RaV2v0o]} 


(54) 
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and 


T - CK.  - u'J  - T.  - -2— 
1 A 1 0 l6Ra 


a2v^0  - aRa1/2  C(m  + l)  vJL 

M 


+ fa  - l^aC2  v„  + fa  + l^aC 


'12 


12  - (m  - l)  aRa1/2  Cv£4 

•fc 

+ i»JL iiW#„  1 (56) 

2 15  J 


where  the  primes  on  the  u-functions  denote  differentiations  with  respect 
to  y and  on  the  v-functions  denote  differentiations  with  respect  to  tj, 
and  the  explicit  forms  of  the  second  derivatives  appearing  in  equations 
(54)  and  (55)  are  presented  in  equations  (B7 ) to  (B13)  of  appendix  B. 

Thus,  solutions  of  the  original  boundary- value  problem  in  terms  of 
u and  T are  known  to  zeroth-  and  first-order  approximations,  that  is, 
neglecting  frictional  heating  and  including  its  effects  to  a first  order, 
respectively.  These  solutions  are,  however,  given  in  terms  of  the  param- 
eter C (recall  that  m is  also  a function  of  C through  eq.  (25)). 
Therefore,  to  relate  C to  the  physical  problem  the  dimensionless  mass 
flow  in  the  channel  is  defined  as 


I 


t 


* 


M 


u0dy 


(56) 


Neglecting  aerodynamic  heating,  equation  (56)  becomes  (using  the  zeroth- 
order  approximation  as  given  by  eq.  (52)) 


a _ a /sin  2r  + sinh  2r\ 
Ra  Ral/4  \posh  2r  + cos  2r/ 


or 


C - 


+ Ra1/*  C(m  +..1)  (l 

V2  V 


sin  2r  - sinh  2r 
cosh  2r  + cos  2r 


)] 


(57) 


cosh  2r  + cos  2r 


2Ad 


a/2  Ra 


(c 

174  in  2r  - sinh  2r 


MRa  a >J2  a/sin  2r  + sinh  2A 
'/  Ka  " Ral/4\cosh  2r  + cos  2tJ 


(58) 
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(For  large  Ra  the  expressions  in  parentheses  in  equation  (58)  reduce 
to  1 and  -1,  respectively. ) Substitution  of  C as  determined  into  equa- 
tions (52)  to  (55)  then  yields  the  velocity  and  temperature  distributions 
for  any  case  if  M,  Ad,  Ra,  a,  and  KA  are  known.  Note,  then,  that 

the  mass  flow  is  an  independent  parameter  of  the  problem,  and  hence  could 
be  due  in  part  to  a forced  flow.  For  a given  configuration  and  fluid  the 
last  four  of  these  parameters  are  specified,  so  that  equation  (58)  re- 
lates C to  the  mass  flow  M,  that  is,  the  temperature  level  and  the 
mass  flow  through  the  channel  are  related,  as  is  only  reasonable  in  flows 
of  this  type. 


Asymptotic  solutions.  - The  solutions  presented  in  the  previous  sec- 
tions  are  valid  (to  the  proper  order  of  approximation)  for  all  values  of 
the  parameters  of  the  boundary-value  problem.  However,  from  physical 
considerations,  it  can  be  seen  that  in  many  practical  occurrences  of  the 
phenomenon  under  consideration,  the  parameter  Ra  may  become  very  large 
(of  the  order  of  104  and  higher) . It  is  therefore  appropriate  to  examine 
the  asymptotic  character  of  the  boundary-value  problem.  To  this  end  it 
is  convenient  to  write  equations  (22)  to  (25)  as 


1_ 

Ra 


_Lv 

u 


KAC  2 

ST 


a 

CRa 


0 


(59) 


u(0)  ■ u(l)  = 0 (60) 

u"(0)  a 1 u”(l)  = m (61) 

where 

u =*  CK^u  (62) 

For  very  large  Ra,  equation  (59)  is  of  the  boundary  layer  type  (see  ref. 
10 ).  Therefore  the  velocity  and  temperature  profiles  will  have  very  large 
gradients  near  the  walls,  and  thus  the  asymptotic  solutions  will  yield  the 
velocity  and  thermal  boundary  lasers  with  essentially  constant  conditions 
given  by  the  inviscid  solution  u^  = a/CRa  in  the  center  of  the  channel 

associated  with  large  Ra  flows.  Hence,  expanding  the  coordinate  normal 
to  the  wall,  as  is  done  in  boundary  layer  theory,  requires  that 

y = RaV^y  (63) 


u » Ra-1/^ 


(64) 
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The  asymptotic  forms  of  equations  (59)  to  (61 ) are 


vyyyy 

C^a  , . 2 a 

+ v - — (V  - —772  - 0 

Ra  y CRa1' 

(65) 

V(0)  = 0 

(66) 

T57(0)-5-{i 

for  the  wall  at  temperature  TWq 
for  the  other  wall  at  temperature  T„ 

1 

(67) 

V(»)  = V^(oo)  - 0 

(68) 

where  the  subscripts  denote  differentiation.  Note  that  the  conditions 
expressed  by  equation  (68)  replace  the  boundary  conditions  at  the  second 
wall  and  require  that  influences  of  one  wall  do  not  affect  the  other. 

Thus  y can  be  considered  as  the  coordinate  normal  to  the  first  (or  left- 
hand)  wall  (that  is,  the  one  corresponding  to  Y = 0)  and  -y  will  be 
the  coordinate  normal  to  the  other  wall,  which  is  at  Y - d.  Hence,  use 
should  be  made  of  the  proper  part  of  equation  (67)  in  each  solution. 

From  equation  (65)  it  can  be  seen  that  the  frictional  heating  effects 
will  be  negligible  for  large  Ra  unless  CK^/Ra  (or  CfPfxd/cp)  is  at 

least  of  unit  order  of  magnitude.  It  should  perhaps  be  pointed  out  here, 
in  contradistinction  to  the  qualitative  discussions  in  references  1 and 
12,  that  the  frictional  heating  is  important  only  if  essentially  the 
ratio  of  K (based  on  any  appropriate  temperature)  to  Ra  is  of  unit 
order  of  magnitude  or  larger,  as  can,  in  fact,  be  verified  in  general. 

Thus  it  should  be  noted  that  the  discussions  in  those  references  hold 
specifically  only  if  Ra  is  of  unit  order  or  smaller  or  if  Ra  does  not 
appear  as  an  explicit  parameter  of  the  problem  (as  in  ref.  1,  for  example). 
For  the  range  of  conditions  and  physical  properties  of  fluids  being  con- 
sidered, it  is  unlikely  that  CK^/Ra  will  be  of  unit  order  for  a flow 

generated  in  a gravitational  field  alone  with  large  Ra.  Therefore,  un- 
less the  natural-convection  flow  is  being  generated  by  a body  force  con- 
siderably stronger  than  gravity,  the  ratio  of  the  volumetric  expansion 
coefficient  to  the  specific  heat  at  constant  pressure  is  unusually  large, 
or  there  is  considerable  forced  flow  (to  increase  C),  the  frictional 
heating  effects  will  not  be  important  for  large  Ra.  (Of  course,  the 
possibility  always  exists  that  some  unusual  fluid  will  be  employed  whose 
physical  properties  are  such  that  CK^/Ra  will  be  of  such  a magnitude 

that  the  frictional  heating  effects  will  be  important  for  large  Ra  even 
in  a gravitational  field.  Liquids  near  their  critical  state  may  be  rep- 
resentative in  this  respect  (see  ref.  13)).  Furthermore,  it  can  be  seen 


<< 
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from  the  same  equation  that  the  effects  of  the  heat  sources  will  he  im- 
portant only  if  a/CRa^/2  is  of  unit  order  of  magnitude  for  large  Ra. 
This  is  physically  reasonable. 

Since  the  frictional  heating  effects  are  negligible  for  large  Ra, 
unless  CK^  is  very  large  these  effects  will  be  neglected  in  this  sec- 
tion. . (A  method  of  successive  approximations  similar  to  that  described 
in  the  previous  sections  could  be  applied  without  difficulty  to  eq.  (65) 
if  these  effects  are  of  consequence.)  Therefore,  letting 

V =*  V - a/CRa1/2,  equations  (65)  to  (68)  become,  respectively, 

V + V = 0 (69] 

yyyy 


Vy£(0)  - m 

V(«)  = v—  - 0 

The  solution  satisfying  equations  (69)  to  (72)  is 


V = -e  ^ hn  sin  + — ~~T~  cos  ) (7 

\ V2  CRa1'2  VV 

To  find  the  temperature  distribution  the  second  derivative  of  equation 
(73)  is  necessary,  and  this  is  given  by 


^"e"y/V5(^slnvi'Sc0£^)  (?4) 

For  large  values  of  the  Rayleigh  number  Ra,  then,  the  dimensionless 
velocity  and  temperature  distributions  are  given  by 


“ “I/i  (v 


(7  + ^) 
\ CRa  ' / 


T a - CKA  (1  - V— ) 
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Solution  of  Special  Case  Simulating  an  Enclosed  Channel 

There  is  considerable  interest  in  the  natural-convection  flow  in  a 
completely  enclosed  rectangular  region.  In  reference  7,  the  natural- 
convection  flow  of  fluids  containing  heat  sources  between  two  parallel 
planes  is  considered  (as  is  the  case  here  also);  but  there,  in  order  to 
simulate  an  enclosed  rectangular  region,  it  is  further  specified  that  the 
net  mass  flow  be  zero  and  that  the  walls  be  at  the  same  temperature.  This 
problem  is  treated  in  an  approximate  manner  in  reference  7 in  that  the 
velocity  distribution  was  postulated  without  regard  to  the  equations  of 
motion  and,  hence,  it  would  be  desirable  for  comparison  purposes  to  obtain 
a more  exact  solution  for  this  special  case  from  the  solutions  found 
herein. 


For  zero  net  mass  flow  in  the  channel  the  parameter  C can  be  deter- 
mined directly  from  the  other  parameters  of  the  problem  and  equation  (58) . 
Thus 


2Ad 


\[2  Ra1 


cosh  2r  + cos  2r\ 


l/4  \sin  2r  - sinh  2r/ 


i)  k 


sin  2r  + sinh  2r 
.cosh  2r  + cos  2r 


H 


(77a) 


where  the  superscript  denotes  the  zero  net  mass  flow  and  the  subscript 
designates  the  order  of  the  approximation.  It  is  interesting  to  note 
from  this  equation  that  for  no  internal  heat  sources  (a  » 0),  zero  net 
mass  flow  in  the  channel  can  be  obtained  when  C ■ 0W  /2Ad.  However,  to 


obtain  the  solution  for  the  special  case  simulating  flow  with  heat  sources 
in  a completely  enclosed  region  with  walls  of  equal  temperatures  from  the 
solutions  presented  in  the  previous  sections,  0W  must  be  zero  in  equa- 
tion (77a),  and  therefore 


H 

8 


,0 

'0 


_ a 

/cosh  2r  + cos 

2r) 

^/5 

/sin  2r  + sinh 

2r\  . x 

V2  Ra1/4 

\sin  2r  - sinh 

2r/ 

_ Ra1/4 

Icosh  2r  + cos 

2r/ 

(77b) 


The  velocity  and  temperature  distributions  for  this  special  case  are  then 


obtained  by  replacing  C in  equations  (52)  and  (54)  by 


as  given  by 


equation  (77b).  Further  physical  significance  of  the  zero  net  mass  flow 
case  can  be  inferred  from  equation  (13).  Integration  of  this  equation 
over  the  channel  cross  section  shows  that  for  no  net  mass  flow  all  the 
heat  generated  internally  (by  heat  sources  and  by  aerodynamic  heating  if 
the  latter  is  significant)  in  a given  cross  section  must  be  transferred 
to  the  walls. 
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Computations  made  for  this  case  including  frictional  heating  yielded 
no  appreciable  deviations  from  the  zeroth-order  results. 


Solutions  for  Case  of  Convective  Inversion 

An  interesting  aspect  associated  with  the  natural-convection 
phenomenon  is  that  the  Grashof  number  can  change  sign;  this  implies  a 
reversal  in  the  flow  direction  and  is  referred  to  as  convective  inver- 
sion (see  p.  109,  ref.  8).  The  sign  of  the  modified  Grashof  number  in 
the  present  paper  can  be  changed  in  one  of  three  ways:  (l)  by  a change 

of  the  sign  of  the  longitudinal  temperature  gradient  dT  /dX  » A,  (2)  by 
a change  in  the  direction  of  the  generating  body  force,  and  (3)  by  a 
change  in  the  sign  of  the  volumetric  expansion  coefficient  as  occurs 
near  the  critical  state  of  a liquid  (see  ref.  13). 

Since  the  modified  Grashof  number  appears  in  the  parameters 

and  Ra  connected  with  the  problem  considered  herein,  the  effects  of 
sign  changes  of  the  Grashof  number  in  the  solutions^should  be  studied. 
From  its  definition,  is  proportional  to  A,  fy  , and  fr  so  that 

only  the  first  will  alter  its  sign.  Note  further,  however,  that  the 
modified  Rayleigh  number  Ra  is  essentially  the  product  of  the  Prandtl 
number  and  a modified  Grashof  number,  that  is,  a Grashof  number  which 
depends  on  the  product  (Ad)  of  the  longitudinal  temperature  gradient  and 
the  distance  between  the  plates.  Hence,  any  one  of  (l),  (2),  or  (3) 
given  in  the  preceding  paragraph  will  lead  to  a change  in  the  sign  of 
Ra.  For  negative  Ra  the  solutions  as  given  in  the  previous  sections 
do  not  apply,  and  hence  the  foregoing  boundary -value  problem  (eqs.  (22) 
to  (25))  would  have  to  be  solved  with  negative  Ra.  These  solutions 
can  be  readily  obtained,  but  it  is  found  that  with  frictional  heating 
neglected  these  solutions  change  character  with  changes  in  Ra  and  that 
there  exist  critical  negative  values  of  Ra  for  which  the  solutions 
become  meaningless.  In  an  attempt  to  explain  these  unusual  results, 
further  interpretation  of  the  problem  must  be  made.  Reexamination  of 
the  meaning  of  negative  Ra  shows  that  not  only  changes  in  the  body 
force  direction  and  sign  of  the  volumetric  expansion  coefficient  but 
also  a change  in  the  sign  of  the  longitudinal  temperature  gradient  A 
can  lead  to  negative  Ra.  If  the  negative  Ra  is  attributed  to  the 
last  cause,  the  physical  interpretation  of  the  unusual  mathematical 
results  pointed  out  becomes  clearer,  because  a negative  A implies  that 
the  fluid  is  being  heated  from  below  and  this  situation  leads  to  a 
"Rayleigh-type"  instability  of  natural -convection  flows  due  to  the 
"piling  of  heavy  fluid  on  lighter  fluid."  Analogous  interpretations,  of 
course,  also  follow  directly  for  changes  in  the  body  force  direction  and 
in  the  sign  of  the  volumetric  expansion  coefficients.  Natural-convection 
flows  heated  from  below  between  horizontal  plates  have  been  studied  ex- 
perimentally in  some  detail  (see  refs.  8,  14,  and  15,  for  example),  and 
it  was  found  that  the  flow  does  indeed  change  character  (into  cellular 
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motion)  for  certain  critical  values  of  the  Rayleigh  number.  Hence,  it 
is  believed  that  the  critical  values  of  Ra  found  from  the  linearized 
analysis  (in  that  the  aerodynamic  heating  was  neglected)  for  negative 
Ra  may  be  analogous  to  those  observed  in  actual  cases.  However,  since 
this  instability  leads  to  these  additional  complications  and  should  be 
further  investigated  (perhaps  using  the  true  nonlinear  eqs.),  the  case 
of  negative  Ra  will  not  be  treated  further  herein.  It  should,  how- 
ever, be  kept  in  mind  that  if  in  an  actual  case  of  the  configuration 
considered  herein  the  Ra  is  negative,  the  flow  and  heat  transfer  will 
not  be  as  predicted  in  this  paper  but  should  be  expected  to  exhibit  a 
behavior  pertinent  to  the  "unstable-type"  flows. 


RESULTS  AND  DISCUSSION 


Velocity  and  Temperature  Distributions 


The  relations  between  the  actual  and  dimensionless  velocities  and 
temperatures  as  determined  from  the  various  transformations  in  the 
analysis  (see  eqs.  (16))  are 


(78) 


(79) 


where  U and  9 denote  the  actual  and  u and  T,  the  dimensionless 
quantities.  For  a given  fluid,  configuration,  heat-source  intensity, 
and  mass  flow,  the  velocity  and  temperature  distributions  can  be  com- 
puted from  equations  (52)  to  (55)  (for  Ra  > 0);  and  for  zero  net  mass 
flow  and  the  walls  of  the  same  temperature,  by  applying  equation  (77b) 
to  equations  (52)  and  (55).  These  computations  will  be  accurate  within 
the  limits  of  the  method  of  solution;  that  is,  for  moderate  and  small 
values  of  Ra  the  solutions  yield  results  of  reasonable  accuracy  for 
small  CK^,  and  for  large  Ra  the  zeroth-order  approximations  or,  even 

more  simply,  the  asymptotic  solutions  will  give  answers  valid  for  all 
CK^.  The  range  of  applicability  of  the  various  solutions  presented 

herein  will  be  discussed  more  fully  subsequently. 

Because  the  solutions  were  obtained  in  the  convenient  forms  (eqs. 
(52)  and  (54),  for  example)  wherefrom  the  qualitative  effects  of  the 
various  parameters  associated  with  the  problem  can  be  studied,  and  since 


3016 


NACA  TN  3141 


21 


tabular  values  of  the  universal  functions  are  presented  to  facilitate 
computations  for  any  specific  case,  no  extensive  detailed  calculations 
will  be  given  covering  the  entire  range  of  values  taken  on  by  the  param- 
eters. Representative  velocity  and  temperature  profiles  were,  however, 
calculated  for  Ka  ■ 10,  m ■ -1,  1,  and  2,  a ■ 0,  10,  and  100,  and 
Ra  ■ 10,  lO^,  1600,  and  10^.  In  addition,  the  parameter  C was  given 
the  value  -1  in  all  the  computations  except  those  for  the  case  simu- 
lating flow  in  a completely  enclosed  region.  In  this  way  the  relative 
influences  of  the  other  parameters  are  just  as  apparent,  but  the  number 
of  computations  is  greatly  reduced.  The  results  of  these  computations 
are  presented  in  figures  2 to  9.  The  contents  of  each  specific  figure 
(numbers  2 to  9)  are  listed  in  the  following  table: 


KA 

■ 10 

u 0' 

U1 

T0' 

T1 

m 

a 

Ra  = 10 

102 

1600 

10* 

10 

102 

1600 

104 

-1 

0 

2(a) 

3(a) 

4(a) 

5(a) 

6(a) 

7(a) 

8(a) 

9(a) 

10 

2(a) 

3(a) 

4(a) 

5(a) 

6(a) 

7(a) 

8(a) 

9(a) 

100 

2(a) 

3(a) 

4(a) 

5(a) 

6(a) 

7(a) 

8(a) 

9(a) 

1 

0 

2(b) 

3(b) 

a4(b) 

5(b) 

6(b) 

7(b) 

a8(b) 

9(b) 

10 

a>b2(b) 

3(b) 

4(b) 

5(b) 

a'b6(b) 

7(b) 

8(b) 

9(b) 

100 

2(b) 

3(b) 

4(b) 

5(b) 

6(b) 

7(b) 

8(b) 

9(b) 

2 

0 

a'b2(c) 

3(c) 

4(c) 

C5(c) 

S'b6(c) 

7(c) 

8(c) 

C9(c) 

10 

2(c) 

3(c) 

4(c) 

5(c) 

6(c) 

7(c) 

8(c) 

9(c) 

100 

2(c) 

3(c) 

4(c) 

C5(c) 

6(c) 

7(c) 

8(c) 

C9(c) 

aIncludes  results  for  u or  T. 
^Includes  results  for  u\2)  or  r(2). 
cIncludes  results  for  ua  or  Ta. 


For  each  triplet  of  parametric  values  (m,  a,  Ra)  the  profiles  were  com- 
puted with  frictional  heating  neglected  (by  eqs.  (52)  and  (54)  and  de- 
noted by  Uq  and  T q on  the  figs . ) , with  frictional  heating  included 

to  a first  approximation  (by  eqs.  (53)  and  (55)  and  denoted  by  Uj_ 

and  T-jJ,  and  in  several  specific  cases  with  frictional  heating  completely 

accounted  for  (by  numerical  solution  of  eqs.  (22)  to  (25)  using  a Card- 
Programmed  Electronic  Calculator  and  denoted  by  u and  T ) . For 
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Ra  ■ 104  the  asymptotic  solutions  (given  by  eqs.  (75)  and  (76)  and 
denoted  by  ua  and  Ta)  are  also  included  for  ■ 10,  m - 2,  and 
a ■ 0 and  100.  The  asymptotic  solutions  were  computed  from  each  in- 
dividual wall  to  the  channel  center  and  then  faired  in  to  Join  smoothly. 
From  these  computations  any  qualitative  trends  obtained  by  examination 
of  the  solutions  can  be  further  substantiated  and,  in  addition,  some 
definition  of  the  range  in  which  the  frictional  or  aerodynamic  heating 
exerts  a large  influence  can  be  made.  Calculations  were  also  made  for 
the  special  case  simulating  a completely  enclosed  region  in  which  there 
is  no  net  mass  flow  and  the  walls  are  at  the  same  temperature  by  applying 
equation  (77)  to  the  appropriate  solutions.  These  curves  are  given  in 
figures  10  and  11.  The  velocity  and  temperature  profiles  (particularly 
for  Ra  «£  104)  axe  qualitatively  similar  to  those  determined  experi- 
mentally in  reference  16. 

Effect  of  different  wall -temperature  configurations  (m  varying) 
and  heat  sources  (a  varying).  - From  equations  (52)  to  (55)  and  their 
related  universal  functions,  it  can  be  seen,  as  expected,  that  an  in- 
crease in  the  wall  temperature  parameter  m or  an  increase  in  the  heat- 
source  parameter  a results  in  larger  velocities  and  higher  tempera- 
tures. These  trends  together  with  that  of  increasing  net  mass  flow,  as 
represented  by  the  area  under  the  u-curves,  with  m and  a can  be  ob- 
served on  figures  2 to  9.  It  can  also  be  seen  from  figures  6 to  9 that 
if  sufficient  heat  is  generated  by  the  heat  sources,  the  direction  of 
heat  transfer  will  be  changed.  In  agreement  with  the  statements  made  in 
the  section  dealing  with  the  zero  net  mass-flow  case,  note  from  fig- 
ures 2(a),  3(a),  4(a),  and  5(a)  that  if  aerodynamic  heating  is  neglected 
for  m **  -1  (since  C = -l)  and  a ■ 0,  there  is  no  net  mass  flow.  In 
general,  the  velocity  distributions  become  more  symmetrical  with  the 
larger  a (see  figs.  2(a),  3(a),  4(a),  and  5(a),  for  example)  because 
the  heat  added  uniformly  by  the  heat  sources  counteracts  any  asymmetry 
imposed  by  the  wall  thermal  conditions. 

Effect  of  the  modified  Rayleigh  number  (Ra).  - Examination  of 
the  solutions  (eqs.  (52)  to  (55))  shows  that  the  velocities  and  tem- 
peratures decrease  with  increasing  values  of  the  modified  Rayleigh 
number  Ra.  This  trend  can  also  be  seen  by  comparison  of  corresponding 
curves  in  figures  2 to  9,  and  even  by  comparing  with  the  curves  in  refer- 
ence 1 which  are  for  Ra  =■  0.  For  large  Ra  it  can  be  seen  from  fig- 
ures 5(b),  5(c),  9(b),  and  9(c)  that  the  velocity  and  temperature  pro- 
files take  on  a "boundary-layer  form."  Asymptotic  solutions  computed 
for  Ra  ■ 104,  m 10,  m ■ 2,  and  a ■ 0 and  100  are  also  presented 
on  figures  5(c)  and  9(c),  and  these  very  closely  approximate  the  more 
exact  solutions.  Hence,  for  large  Ra  the  asymptotic  solutions  can  be 
employed  to  yield  reasonable  results  much  more  simply.  For  the  case 
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where  m » -1  and  a a 0,  increasing  the  modified  Rayleigh  number  changes 
the  temperature  distribution  from  essentially  the  conduction  profile 
(that  is , an  almost  linear  distribution)  at  Ra  = 10;  the  increased 
effect  of  the  convection  is  then  apparent  for  the  larger  values  of  the 
Rayleigh  number.  It  is  interesting  to  note  that  for  Ka  ■ 10,  a ■ 10, 
and  m ■ -1,  1,  and  2,  changes  in  the  Rayleigh  number  can  so  affect  the 
temperature  distributions  that  the  heat  flow  direction  from  one  or  both 
the  walls  can  be  altered.  (Compare  corresponding  parts  of  figs.  6 
to  9.)  This  point  will  be  more  graphically  portrayed  in  the  subsequent 
discussion  of  Nusselt  numbers. 

The  velocity  and  temperature  distributions  (see  figs.  10  and  ll) 
for  the  special  case  considered  herein  of  zero  net  mass  flow  and  walls 
at  the  same  temperature  are  not  in  general  appreciably  altered  in  shape 
by  increases  in  Ra  although  the  velocity  peaks  vary  inversely  with  the 
Rayleigh  number.  The  shape  of  the  velocity  and  temperature  profiles  is 
seen  to  be  qualitatively  the  same  as  that  assumed  in  reference  7. 

Effect  of  frictional  heating.  - By  comparing  the  profiles  presented 
in  figures  2 to  9 computed  by  neglecting  frictional  heating  (denoted  by 
the  subscript  zero)  with  those  computed  including  the  aerodynamic  heat- 
ing to  a first  order  (denoted  by  the  subscript  unity),  the  effect  of  the 
aerodynamic  heating  on  the  velocities  and  temperatures  can  be  studied. 
Numerical  solutions  obtained  of  the  complete  boundary-value  problem 
(eqs.  (22)  to  (25))  in  which  the  frictional  heating  was  entirely  taken 
into  account  are  also  included  (with  no  subscripts)  (see  preceding  table) 
on  figures  2(b),  2(c),  4(b),  6(b),  6(c),  8(b),  10,  and  11  for  comparisons 
with  the  approximate  solutions. 

In  accord  with  the  discussion  on  the  asymptotic  solutions,  it  can 
be  seen  that  when  KA  is  small  compared  with  Ra  (recall  that  C =»  -1 
in  these  calculations)  the  aerodynamic  heating  effects  are  negligible. 
Since  no  computations  were  made  herein  for  KA  > 10,  the  computations 
made  for  Ra  = 1600  and  10 ^ show  no  deviation  between  the  zeroth  and 
first  approximations,  and  these  are  also  coincident  with  the  numerical 
solution  (see  figs.  4(b)  and  8(b)).  Hence,  in  the  range  KA/Ra«  1, 
the  zeroth-order  approximations  will  yield  accurate  results;  if,  in 
addition,  Ra»  1,  the  asymptotic  solutions  provide  a simple  means  of 
obtaining  the  velocity  and  temperature  profiles.  Note  that  the  param- 
eter KA  serves  merely  as  a scale  factor  in  the  zeroth-order  solutions. 

In  the  range  where  K A and  the  Rayleigh  number  are  of  the  same 

order  of  magnitude,  the  frictional  heating  affects  the  results  to  greater 
or  lesser  degree  depending  on  the  particular  amount  of  heat  addition  as 
specified  essentially  by  the  parameters  m and  a (figs.  2,  3,  6, 
and  7)  and,  hence,  the  first-order  approximations  should  be  employed  in 
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this  range.  Comparison  of  the  first-order  solutions  with  the  several 
numerical  solutions  of  the  entire  boundary-value  problem  shows  close 
agreement  for  the  cases  computed.  Hence,  unless  the  conditions  are 
more  severe  than  the  most  extreme  conditions  for  the  range  of  parameters 
considered  herein  as  represented  by  K^/Ra  ■ 1,  m = 2,  and  a ■ 100, 
the  first-order  solutions  will  yield  results  which  include  the  effect 
of  aerodynamic  heating  of  reasonable  accuracy.  In  this  range  where  its 
effects  are  important,  this  frictional  heat,  of  course,  acts  Just  as  do 
the  heat  sources  and  leads  to  increased  velocities  and  altered  tempera- 
ture profiles  and,  consequently,  different  heat-transfer  rates  (figs.  2, 
3,  6,  and  7).. 

For  the  special  case  simulating  flow  in  a completely  enclosed  re- 
gion, the  aerodynamic  heating  did  not  affect  the  results  appreciably  over 
the  range  of  parameters  under  consideration.  A numerical  solution  com- 
pletely including  the  effects  of  aerodynamic  heating  was  made  for 
Ra  a 10;  it  can  be  seen  in  figures  10  and  11  that  this  solution  coin- 
cides with  the  zeroth-order  solution,  which  neglects  the  effect  of  fric- 
tional heating. 

The  complete  consideration  of  frictional  heating  (as  in  the  numeri- 
cal solutions)  for  the  problem  discussed  herein,  Ju6t  as  for  the  case 
(essentially  Ra  ■ 0)  reported  in  reference  1,  leads  to  the  two  results 
(a)  that  there  exists  a critical  set  of  conditions  beyond  which  no  solu- 
tions exist,  and  (b)  that  where  solutions  exist  there  are  two  solutions 
for  every  set  of  admissible  parametric  values.  Examples  of  these  second 
solutions  are  presented  in  figures  2(b),  2(c),  6(b),  and  6(c)  (denoted 
by  u(2)  and  T(2)),  and  it  can  be  seen  that  the  velocities  are  more 
than  10  times  as  large  as  the  first  solutions  and  the  temperatures  are 
much  greater  than  the  corresponding  first  solution  temperatures.  These 
last  unusual  results  cannot  be  predicted  from  the  solutions  obtained  by 
successive  approximations  as  described  herein,  but  are  found  from  numeri- 
cal solutions  obtained  by  means  of  a Card-Programmed  Electronic  Calcu- 
lator. At  present  the  significance  of  the  second  solutions  is  not  ex- 
plained, although  it  is  felt  that  they  are  intimately  connected  with  the 
unique  regenerative  action  of  the  frictional  heating  in  natural  convec- 
tion. The  existence  of  the  critical  conditions  appears  to  be  similar 
to  the  thermal  choking  phenomenon. 


HEAT  TRANSFER 
Nusselt  Numbers 

The  heat -transfer  coefficients  for  the  natural-convection  phenomenon 
treated  here  can  be  expressed  in  terms  of  Nusselt  numbers.  The  Nusselt 
number  is  here  defined  as 


3016 


where  the  double  subscript  signifies  that  the  temperature  gradient  is 
to  be  evaluated  at  either  Y ■ 0 or  Y - d,  depending  on  the  wall  under 
consideration. 


In  terms  of  the  dimensionless  quantities, 


The  temperature  gradient  can  be  found  from  the  zeroth-order  solutions, 
and  the  Nusselt  numbers  can  be  computed  on  this  basis  from 


L__J  1 - a(sinh  2r  + sin  2r) 

Z-Jz  Ra  ' l cosh  r + cos  r - 1 _ 

(sinh  2r  - sin  2r)  + — = — (sinh  2r  + sin  2r)> 

J 2(cos  r - cosh  r)  J 


2(cos<: 


2r  + sin  2r] 


where  the  first  subscript  denotes  that  zeroth-order  approximation  is 
used,  and  the  second  denotes  the  wall  with  which  the  Nusselt  number  is 
associated.  (Eqs.  (81)  and  (82)  are,  of  course,  specifically  for  Ra 
positive.)  Note  that  these  zeroth-order  Nusselt  numbers  are  independent 
of  K^.  When  CKa  is  of  the  same  order  of  magnitude  as  Ra,  the  zeroth- 

order  approximation  has  been  shown  to  be  inaccurate;  therefore  improved 
Nusselt  numbers  can  be  obtained  by  using  the  appropriate  T (first- 
order  approximation)  solutions  and,  for  the  same  conditions  as  for 
equations  (8l)  and  (82) , the  Nusselt  numbers  can  be  computed  from 

K r 

Nuio  - Nuqq  + -A-  «A\  - CoRa1/2  (m  + ljOg  + 

64 ft  Ra''  L 

+ c2lm  ■t..1?.2Rafl4  + c(m  - lJaRa1/^  . gffcsLz.  1)  (83) 
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«uU  - *01  - -pSyj  - Ca^/Hm  * 1 »2  4 £!iE^5»Q3 

64V2  Ra  ' L- 

♦ c2(‘  + I)2Ran4  + C(b  - l)aB.V%fe  4 C-2fr22--.  11  fcflg]  (84) 

where  the  various  2.^  in  equations  (83)  and  (84)  are  given  by  equations 
(B43)  to  (B48)  in  appendix  B. 

Computations  of  the  Nusselt  numbers  were  made  over  a range  of 
values  of  Ra  from  equations  (81 ) to  (84)  for  = 10,  m = -1,  1,  and 
2,  and  a ■ 0 and  10  where  again  C was  taken  to  be  -1,  and  the  re- 
sults are  presented  in  figures  12  and  13  for  the  wall  at  Y » 0 and 
Y = d,  respectively.  The  figures  show  that  for  the  wall  at  Y ■ 0 the 
Nusselt  numbers  decrease,  in  general,  with  increasing  Ra,  and  for  the 
wall  at  Y * d the  Nusselt  numbers  increase  with  the  modified  Rayleigh 
number  except  for  the  case  where  m = -1.  For  = 10  and  a ■ 10, 
the  Nusselt  number  changes  sign  with  increases  in  Ra  because  of  the 
variation  of  the  temperature  profile  with  Ra,  a6  was  previously  noted. 

Figures  12  and  13  also  demonstrate  clearly  the  effect  of  the  aero- 
dynamic or  frictional  heating.  This  effect,  in  accordance  with  all  that 
preceded,  is  extremely  pronounced  for  low  values  of  the  Rayleigh  number 
(that  is,  when  K/Ra  is  of  unit  order  of  magnitude). 

Flow  in  an  enclosed  region.  - For  the  special  case  simulating  flow 
in  a completely  enclosed  region  (M  - 0,  m = 1,  and  a / 0),  the  calcula- 
tions for  the  temperature  profile  were  extended  over  a larger  range  of 
Ra  and  plotted  in  figure  14  as  the  ratio  $ of  the  temperature  differ- 
ence to  that  for  pure  conduction,  as  was  done  in  reference  7.  The  con- 
duction temperature  difference  used  in  ® is  the  channel  center-to- 
wall  difference  subject  to  uniform  heat  generation  by  sources  and  is 
equal  to  aK^/8.  The  temperature  profiles  computed  in  an  approximate 
manner  in  reference  7 are  compared  with  those  computed  more  exactly  by 
the  method  reported  in  this  paper,  and  it  can  be  seen  that  the  dis- 
crepancy becomes  quite  apparent  for  values  of  Ra  = 104  and  above. 

The  Nj  in  reference  7 is  related  to  Ra  by  Ra  **  — Nj;  hence,  for 

most  liquids  Nj  and  Ra  sure  identical.  If  r / 1,  T is  merely  a 
scale  factor.  It  can  be  seen  from  figure  14  that  for  Ra  < 104,  the 
temperature  gradients  at  the  wall  are  all  almost  identical,  and  hence, 
even  if  the  temperature  profiles  themselves  were  not  identical,  the 
heat  transfer  computed  by  the  two  methods  would  be  in  reasonable  agree - 
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ment.  The  variations  of  the  dimensionless  temperature  variable 
(ratio  of  center-to-wall  temperature  difference  to  that  for  pure  con- 
duction) as  used  in  reference  7 with  Ra  as  given  by  the  two  methods 
are  also  compared  in  figure  15,  and  hence  the  quantitative  limits  of 
the  approximate  method  can  be  seen. 
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CONCLUDING  REMARKS 

An  analysis  was  made  of  the  flow  subject  to  body  forces  between  two 
parallel  plane  surfaces  oriented  in  the  direction  of  the  generating  body 
force  along  which  the  temperature  is  specified  to  vary  linearly.  The 
solutions  for  natural  convection  and  those  for  combined  natural  and 
forced  convection  were  found  to  be  essentially  the  same.  It  was  found 
that  a modified  Rayleigh  number  (product  of  Prandtl  and  modified  Grashof 
numbers)  in  addition  to  the  parameter  Ka  was  of  significance  in  this 
problem.  Solutions  for  the  velocity  and  temperature  distributions  are 
given  in  terms  of  "universal"  tabulated  functions.  Detailed  velocity 
and  temperature  profiles  were  computed  and  it  was  found  that,  in  general, 
the  velocity  and  temperature  differences  increase  with  the  wall  tempera- 
ture parameter  and  with  additional  heat  due  to  heat  sources.  The  veloc- 
ities and  temperatures  decrease  with  increasing  values  of  the  modified 
Rayleigh  number.  When  the  ratio  of  CKa  to  the  modified  Rayleigh  num- 
ber is  of  unit  order  of  magnitude,  the  frictional  or  aerodynamic  heating 
appreciably  affecns  the  velocity  and  temperature  distributions.  Asymptotic 
solutions  for  large  Ra  are  presented  which  make  computations  in  this 
range  relatively  simple.  For  any  given  set  of  the  parameters,  complete 
consideration  of  frictional  heating  implied  the  existence  of  two  flow  and 
heat-transfer  states  and  implied  that  no  solution  exists  beyond  certain 
critical  values. 

Consideration  was  given  to  a special  case  simulating  the  natural- 
convection  flow  of  fluids  with  heat  sources  in  a completely  enclosed 
region  with  the  walls  at  the  same  temperature.  Computations  from  the 
solutions  for  this  special  case  demonstrated  that  an  approximate  method 
developed  in  another  paper  should  yield  reasonably  accurate  results  as 
long  as  the  modified  Rayleigh  number  is  less  than  104.  The  effects  of 
aerodynamic  heating  were  found  to  be  negligible  for  this  case. 

Study  of  the  convective  inversion  aspect  of  the  present  problem  led 
to  the  inference  that  for  negative  values  of  the  Rayleigh  number,  addi- 
tional complications  arise  because  of  an  instability  (due  to  heating  from 
below)  of  the  flow  which  must  be  more  thoroughly  investigated. 


Lewis  Flight  Propulsion  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Cleveland,  Ohio,  December  29,  1953 
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APPENDIX  A 
SYMBOLS 

The  following  notation  is  used  in  this  report: 

A longitudinal  temperature  gradient 

n 

a^  constants  in  eq.  (39);  i = 1,  2,  3,  4 

-*V 

constants  defined  "by  eqs.  (B14)  to  (B35) ; i = 1,  2,  3,... 22 

C constant  in  eq.  (17) 

constants  defined  by  eqs.  (B36) ; i » 1,  2,  3,  ...6 

C®  constant  defined  by  eq.  (77b) 

C constant  in  eq.  (14b) 

Cp  specific  heat  at  constant  pressure 

cv  specific  heat  at  constant  volume 

Dj_  constants  defined  by  eqs.  (B37);  i m 1,  2,  3,  ...10 

d characteristic  length  (specifically  distance  between  plates) 

E constant  defined  by  eq.  (50) 

constants  defined  by  eqs.  (B38) ; * _ 1,  2,  3,  ...6 

e constant,  (i+l) 

F constant  defined  by  eq.  (5l) 

f^  negative  of  X-component  of  body  force  per  unit  mass 
G Green's  function 

4/2 

Gr^  modified  Grashof  number,  pf^  Ad  /y 
g constant,  (i-l) 


}'  •«**'.  3, 
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h heat-tranefer  coefficient 

hi  constants  defined  by  eq.  (B39) ; i = 1,  2,  3,  4 
J constant,  Ra^/^CK^/64 

Pfxd 

dimensionless  parameter,  Pr  Gr^  “ — 

k thermal-conductivity  coefficient 

^ constants  defined  by  eqs.  (B40);  i = 1,  2 

M dimensionless  mass  flow 

Mi  constants  defined  by  eqs.  (B4l);  i = 1,  2,  3,  4 
m constant  defined  by  eq.  (25) 

m constant  defined  by  eq.  (64) 

Nu  Nusselt  number,  hd/k 

Nj  modified  Rayleigh  criterion  as  given  in  ref.  6,  yRa 
P pressure 

Pr  Prandtl  number,  cp|i/k 

Q heat  due  to  heat  sources 

Qi  constants  defined  by  eqs.  (B42);  i = 1,  2,  3,... 8 
R constant,  Ra/l6 

Ra  modified  Rayleigh  number,  i Pr  Gr^ 

r constant, 

T,  T*  temperature 
U velocity 

u dimensionless  velocity 

u dimensionless  velocity,  u/CK^ 


■ - 

„ rj. 
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dimensionless  velocity  defined  "by  eq.  (6l) 


dimensionless  velocity,  V 


dimensionless  velocity,  v/64 
dimensionless  velocity  defined  by  eq.  (26) 
longitudinal  coordinate 
transverse  coordinate 


dimensionless  transverse  coordinate 


dimensionless  transverse  coordinate 


dimensionless  heat-source  parameter,  Qd/kA 


coefficient  of  volumetric  expansion,  p 


ratio  of  specific  heats 
dimensionless  coordinate,  2y  - 1 
temperature  difference,  T " Two 

constant , -a KpJ 64 

absolute  viscosity  coefficient 
kinematic  viscosity  coefficient 
dummy  variable 
density 

dimensionless  temperature  difference 
dimensionless  temperature  difference,  8t/oK^ 

dimensionless  center-to-vall  temperature  difference,  8('T)y6-i/^aK/v 
constants  in  eqs.  (80)  and  (8l)j  i » 1,  2,  3, ...6 


3S.V. 
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Subscripts: 

a asymptotic  solution 

ca  complementary  antisymmetric  solution 

cs  complementary  symmetric  solution 

I inviscid  solution 

n order  of  approximation 

P particular  solution 

wq  conditions  at  y = 0 

w^  conditions  at  y = 1 

0 zeroth-order  approximation 

1 first-order  approximation 
Superscript : 

(2)  second  flow  and  heat-transfer  state 


t 
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APPENDIX  B 


FIRST-CRDER  SOLUTIONS 


TCie  explicit  forms  of  the  first-order  solutions  to  be  used  in  equa- 
tion (46)  ares 

v10  = — bfBl  cosh  rTl  cos  rr)  + B2  sinh  rt)  sin  rrj  + B3  cosh  2rrj  cos  2rrj 
Ra  ' u 

+ B4  sinh  2rr\  sin  2rT]  + B5(cos  2rr\  - cosh  2rii)  + 15  B^J  (Bl) 

V11  = — 572  O®6  cosh  rT)  cos  rT)  + B7  sinh  rTl  sln  rTl  " b4  cosh  2rri  cos  2ri] 
R&  / 

+ B3  sinh  2rT)  sin  2rij  - 15  (B2) 

v12  ■ Ra3/2[B8  cosh  rrj  cos  r»i  + Bg  sinh  rTi  sin  rrj 
+ B10  cosh  2 rr)  cos  2 rr)  + B1:L  sinh  2 rrj  sin  2rr) 

+ B12(cosh  2rrj  + cos  2rrj)  - 15  B1Q]  (B3) 

v13  = ,a3/2[B13  cosh  rTl  cos  rn  + B14  sinh  rrj  sin  rrj 
- Bj  cosh  2rrj  cos  2rr\  - B4  sinh  2tt\  sin  2rr) 

+ B5  (cos  2rrj  - cosh  2rrj)  - 15  B^J  (B4) 

*14  = 3/2^15  s^-n^1  rTl  cos  rr)  + B^g  cosh  ri)  sin  ri] 

+ B17  sinh  2rrj  cos  2rrj  + B10  cosh  2tt\  sin  2rrj 

+ Big  sin  2rn  + B2Q  sinh  2rij|  (35) 


> v 


,J-  o 

■■  . - ' ‘ 
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V15  = p 3/2 21  sin^  rT*  C0S  rT^  + ®22  COSh  rT'  8*n  rT* 
- Bj_q  sinh  2rq  cos  2rr)  + B]_y  cosh  2rq  sin  2rr) 


+ B20  sin  2rq  - B^g  sinh  2rqJ  (B6) 

The  constants  B^(i  = 1,  2,  ...  22)  appearing  in  the  preceding  equations 

are  readily  computed  for  a given  Fa.  In  the  subsequent  section  of  this 
appendix  these  constants  are  written  explicitly  in  a form  suitable  for 
reasonably  rapid  computation. 

The  second  derivatives  appearing  in  equations  (54)  and  (55)  are 
ti  -RaV^  r 

v"  = = 5 — Icosh  r sin  r cosh  rq  sin  rq 

4(cos^r  - cosher) 

+ sinh  r cos  r sinh  rq  cos  rqj  (B7) 

v10  = Sa[B2  C0Sk  rT'  C0S  rT*  " B1  Sinl1  rT)  sin  rT^ 

+ 4B4  cosh  2rrj  cos  2rq  - 4Bj  sinh  2rq  sin  2rq 

- 2B5(cos  2rq  + cosh  2rq)J  (B8) 

V11  = Sa[B7  cosh  rq  cos  rq  - Bg  sinh  rq  sin  rq 

+ 4B3  cosh  2rq  cos  2rq  + 4B4  inh  2rq  sin  2rqJ  (B9) 

v12  = 4Ra[^9  cos^  rT)  cos  rTl  ” ®8  rTl  sin  rTl 

+ 4Bj2_  cosh  2rq  cos  2rq  - 4B-j_q  sinh  2rq  sin  2rq 
+ 2B^g(cosh  2rq  - cos  2rq)j  (BIO) 

V13  = 4^ai?l4  COsl1  rT^  C0S  rT^  “ ®13  s^nh  r7l  sin  rTl 

- 4B4  cosh  2rq  cos  2rq  + 4B3  sinh  2rq  sin  2rq 

- 2Bg(cos  2rq  + cosh  2rq)]  (Bll) 


fH  : 

Pi  : 


34 


NACA  TN  3141 


v£4  = 4^["B15  cosh  rT3  sin  rTl  + B16  sinil  rT)  008  rTl 

- 4B yj  cosh  2rrj  sin  2rrj  + 4B^q  sinh  2rrj  cos  2rrj 

- 2B^g  sin  2rrj  + 2B2o  sinh  2rrjJ  (B12) 

v15  = 4^^"®21  cos^  rr)  sin  rTl  + ®22  rT)  cos  rTJ 

+ 4Bi8  cosh  2rrj  sin  2rr)  + 4Bi7  sinh  2rr)  cos  2rrj 

- 2B2q  s*n  2r,l  “ ^B19  s*11*1  2rrjj  (B13) 

Constants  for  First-Order  Approximation 

The  constants  appearing  in  equations  (Bl)  to  (B13)  are  written  in 
an  expeditious  form  for  computing  as  follows: 


B1  = ‘ i^[("3Dl+D2)El  + (D3-3D4)E2  + 2(D3+D4)E3  + 2(Di+D2)E4 


+ (D3+D4)E5  - (Dj^+DgjEg 

512  hlh2 
15  Mx 

(B14) 

- - ±[(D3-3D4)E1  + (3D1-D2)E2  + 2(D1+D2)E3  - 2(D3+D4)E4 

- (l^+DgjEg  - (D3+D4)Egj 

128  (2hl2  + h22^ 

r 15  mx 

(B15) 

32  hlh2 
= ‘ 15  Mx 

(B16) 

16  ^hl2  " h2^ 
- T5 

(B17) 

16  (hl2  + h22) 
5 M, 

(bis) 

B6  * + <V3D4>E1  + <V°4)E6  + <WE5 

• rlf  ■ hi2)]  (B19) 


<> 

>:> 


\ 
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(j' 


>v- 


1 

3 

3 


B7  = -D3+3D4 )Eg  + (3D1-D2)E1  - (D1+D2)E6  + (D3+D4)E^j- 


128  hjhg 


15Mi 


(B20) 


B8  = ■ 4[(VV2D4)C1  + (2D1+D7-Dq)C2  - 2(D5+D6)C3  - 2(D7+Dq)C4 

, 448  h h,  32(Mc-fM.) 

+ (D5-D6+2D4)C5  + (-2D1+D7-D8)C6J  + (B21) 

[(2D1+D7-Dq)C1  + (-D5+D6+2D4)C2  - 2(D7+D8)C3  + 2(D5+Dg)C4 


B9=  -Z 


/ % - 64(Mr-M4) 

+ (-2D1+D7-De)C5  + (-D5+B6-2D4)C6]+  4 


32  h-^2 
Bl°  = "l5M^ 


(B22) 

(B23) 


Bn  = . 


16(M5-M4) 


B 


12  = 


ISMg 
16(Mg+M4 ) 


5Mo 


(B24) 

(B25) 


B13  " L1t'3Dl+D2)El  + (D3-3D4)E2  - 2(D3+D4)E3  - E^+D^E^ 
/ \ / -1  512  h-iho 

+ (d3+d4)e5  - (d1+d2)e6J- 


15M, 


(B26) 


Bl4  = ^tD3-3D4)El  + (3D1-D2)E2  - 2(D1+D2)E3  + 2(D3+D4)E4  - 

- (D3+D4)Eg]  - ^ [hj2  + 2hg2]  (B27) 


B15  " ^[(2Dl+D7-D8)Ql  + (D5-D6-2D4)Q2  " (D5+De)Q3  + (d7+D8)Q4 

" (°9  sinh  2r)Q5  " (D10  slnh  2r)<*6  - (D9  sln  2r  )Q7  - (D10  sin  2r] 

+ i5M3[-20h3  + 12h4  + 18  sinh  2r  - 6 sin  2r]  (B28) 


. Q 


>- 


■ ¥$  , 


v 
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B16  = ^•[(-D5+I)6+2D4)<5l  + (2Di+C7-D8)Q2  + (D7+Dq)Q3  + (Dg+Dg)^ 
- (D10  sinh  2r)Q5  + (Dg  sinh  2r)Qg  - (D10  sin  2r)Q? 

+ (Dg  sin  2r)QgJ+  Yj|^[l2h3  " 20h4  ” 6 sinh  2r  + 18  sin  2rJ 


Bl7  ~ 15M, 


318  ~ 15K 


Bn  Q = 


8 sinh  2r 


320  = “ 


8 sin  2r 


B21  * L2[(-2Di-D7+D8)Q2  + (D5-D6-2D4)Q1  + (Dg+Dg)^  + (Dy+DgjQj 
+ (Dg  sinh  2r)d7  + (D10  sinh  2r)Qe  - (Dg  sin  2r)Q5 

- (D10  sin  2r)QgJ+  ~“[20h4  + 12h3  + 18  sin  2r  + 6 sinh  2rJ 

I 

B22  = L^[(I)5'D6_2D4)Q2  + (2D1+D7"D8)Q1  “ (VVQ4  + 

+ (D10  sinh  2r)Q?  - (Dg  sinh  2r)Q0  - (d10  sin  2r)Q5 

+ (Dg  Sin  2r)QgJ-  jgjj^[l2h4  + 20h3  + 6 sin  2r  + 18  sinh  2rJ 


jfr  ■ 

51  jh  , 

I : : 
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where 


2 4 

= - — sinh  2r  sin  2r  - ^ sinh  r sin  r 


2 4 

C2  = j cosh  2r  cos  2r  - 2 + - cosh  r cos  r 


2 

5 


C3  = - — cos  2r  - — cosh  2r  + ■?  cosh  r cos  r 


4 

5 


V (B36) 


4 4 4 

C4  - - — cos  2r  + — cosh  2r  + — sinh  r sin  r 


2 12 
C5  = - sinh  2r  sin  2r  + sinh  r sin  r 


12 


Cg  = — cosh  2r  cos  2r  + 2 - -g-  cosh  r cos  r 


Dj_  = cosh2r  cos2r  sinh  r sin  r 


Do  = sinh°r  sin°r 


D*  = cosh°r  cos°r 


D4  = sinh2r  sin2r  cosh  r cos  r 


Dq  = cosh  r cos^r  sinh2r 


D6  = cosher  cos  r sin2r 


D7  = sinh^r  sin  r cos2r 


Do  = cosh2r  sin^r  sinh  r 


Dg  = sinh  r cos  r 


Dio  = cos^  r sin  r 


(B37) 


{ i. 


H j 


X/ 


' **  l 
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2 8 
Eq_  = — cosh  2r  cos  2r  + 2 - ^ cosh  r cos  r 


2 8 
E2  = g-  slnh  2r  sin  2r  - -g  sinh  r sin  r 


2 2 Q 

E*  = - — cos  2r  + - cosh  2r  - - sinh  r sin  r 

° o o o 


4 4 8 

E 4.  - - t cos  2r  - r-  cosh  2r  + — cosh  r cos  r 

* o o b 


2 Q 

E5  = - g sinh  2r  sin  2r  + — sinh  r sin  r 


2 q 

E0  = - - cosh  2r  cos  2r  + 2 - - cosh  r cos  r 


h^  = cosh  r cos  r 


h2  = sinh  r sin  r 


h3  = cosh  r sinh  r(l  - 2 cos2r) 
h4  = - sin  r cos  r(2  cosh2r  - l) 


L-l  = (cosh2r  + cos2r  - l)^ 

L2  = (cosh2r  - cos2r)2  (cosh2r  + cos2r  - l) 


M-|_  = (cosh  r + cos^r  - l)2 
M2  = (cosh2r  - cos2r)2 

M3  = (cosh2r  + cos2r  - l)  (cosh2r  - cos2r) 
M A = cosh2r  sin2r 


Me  = sin^r  cos^r 


¥t 
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t ' 

f 


2 4 

Ql  = ^ cosh  2r  sin  2r  - ^ cosh  r sin  r 

2 4 

Q2  = j s^r^1  cos  2r  - — sinh  r cos  r 

2 4 

Q3  = — sinh  2r  cos  2r  - — cosh  r sin  r 


2 4 

Q4  = - -r  cosh  2r  sin  2r  + - sinh  r cos  r 


Q5  = - ~ sin  2r  + ^ sinh  r cos  r - ^ cosh  r sin  r 


Q5  = - — sin  2r  + — cosh  r sin  r + — sinh  r cos  r 


2 6 2 

Q7  = - — sinh  2r  + — cosh  r sin  r - — sinh  r cos  r 

• D D D 


4 6 2 

Qg  = — sinh  2r  - - sinh  r cos  r - -g  cosh  r sin  r 


> (B42) 


Nusselt  Number  Constants 

The  constants  appearing  in  equations  (83)  to  (84)  are  given 
explicitly  as 

Qa.  = (Bg-Bi  )sinh  r cos  r - (B^+B2)cosh  r sin  r + 8(B4-B3)sinh  2r  cos  2r 

- 8(B3+B4)cosh  2r  sin  2r  - 4B3(sinh  2r  - sin  2r)  (B43) 

Q2  = (B7-B  g)sinh  r cos  r - (B6+B7)cosh  r sin  r + 8(B3+B4)sinh  2r  cos  2r 


+ 8(B4-B3)cosh  2r  sin  2r 


(B44) 


"3  = (Bg-Bg)sinh  r cos  r - (Bg+Bg)cosh  r sin  r + afB^-B^Jsinh  2r  cos  2r 
- 8(B^0+D11)cosh  2r  sin  2r  + 4Bl2(sinh  2r  + sin  2r)  (B45) 


\ i 


, ■ ■ W 

- yj: 


i<5«Ns-'^4 


V; 
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^4  = (B14"B  13)sinh  r cos  r - (B13+B^4)cosh  r sin  r 
+ 8(B3-B4)sinh  2r  cos  2r  + 8(B3+B4)cosh  2r  sin  2r 

- 4Bg(sinh  2r  - sin  2r)  (B46) 

^5  = (B16"B  3_5)cosh  r cos  r - (B^5+B^g)sinh  r sin  r 

+ 8(B]_8-B1’7)cosh  2r  cos  2r  - 8(Bi7+Bpg)sinh  2r  sin  2r 

- 4B^_g  cos  2r  + 4B20  cosh  2r  (B47) 

= (B22-B2x)cosh  r cos  r - (B2}+B22)sinh  r sin  r 

+ 8(B1-7+B2_g)cosh  2r  cos  2r  - 8(B^-Bi8)sinh  2r  sin  2r 

- 4B2o  cos  2r  - 4B-^q  cosh  2r  (B48) 
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TABLE  II.  - FIRST  ORDER  UNIVERSAL  FUNCTIONS 
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(a)  m - -1. 

Figure  2.  - Dines* looless  velocity  distributions  for  various 
heat-source  peraneters  vlth  Ra  ■ 10,  KA  ■ 10,  C « -1. 
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Figure  2.  - Continued.  Dimensionless  velocity  distributions 
for  various  heat-source  parameters  with  Ra  * 10,  Ka  =>  10, 
= -1 
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y - Y/d 
(c)  m ■ 2. 

Figure  2.  - Concluded.  Dimensionless  velocity  distributions 
for  various  heat-source  parameters  with  Ra  * 10,  K*  ■ 10, 
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MCA  TN  3141 


(a)  m » -1. 


Figure  3.  - Dimensionless  velocity  distributions 
heat-source  parameters  with  Ra  « 10  , KA  « 10, 


for  various 
C - -1. 


Dimensionless  velocity 


NACA  TN  3141 


49 


Figure  3.  - Continued.  Dimensionless  velocity  distributions 


for  various  heat-source  parameters  with  Ra 


.6 


y » T/a 

(c)  N - 2. 


Figure  3.  • Concluded.  Dlaenalocleaa  velocity  dlatrlbutlooa 
for  verloua  heat-aource  peraaetera  with  B*  ■ 102,  K.  « 10 
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Figure  4.  - Dimensionless  velocity  distributions  for  various 
heat-source  parameters  with  Ha  * 1600,  - 10,  C = -1. 
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Figure  5.  - Continued.  Dimensionless  velocity  distributions 
for  various  heat-source  parameters  with  Ra  ■ 10,  K.  - 10, 
C - -1.  * 
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(c)  m ■ 2. 

Figure  5.  - Concluded.  Dimensionless  velocity  distributions 
for  various  heat-source  parameters  with  Ra  - 10  , K.  - 10, 
C » -1.  A 


3016 


(t>)  m = 1. 


Figure  6 . - Continued.  Dimensionless  temperature  distribu 
tions  for  various  heat-source  parameters  with  Ra  = 10, 

- 10,  C = -1 
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Figure  7.  - Continued.  Dimensionless  temperature  distribu- 
tions for  various  heat-source  parameters  with  Ra  = 10  , 
Ka  - 10,  C - -1. 
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y = Y/d 


(a)  m ■ -1. 

Figure  8.  - Dimensionless  temperature  distributions  for  vari- 
ous heat-source  parameters  with  Ra  = 1600,  ■ 10,  C = -1. 


(c)  m ft  2. 


Figure  8.  - Concluded.  Dimensionless  temperature  distribu 

tions  for  various  heat-source  parameters  with  Ra.  * 1600 
K.  - 10,  C » -1. 
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Figure  9.  - Continued.  Dimensionless  temperature  distribu- 
tions for  various  heat-source  parameters  with  Ra  ■ 10  , 

Ka  - 10,  C - -1. 
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Figure  10.  - Dlaenelonless  velocities  for  special  case  of 
M ■ 0,  a • 1,  KA  ■ 10,  sod  a ■ 10.  (Hote:  For  R*  a lo 
the  nuaerlcal  solution  u coincides  with  given  curve.) 
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